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Abstract 

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization 
due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, 
driver of biodiversity loss. Our study system, the Alpine whitefish {Coregonus spp.) species complex is a classic example 
of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of 
similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by 
anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the 
most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this 
radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent 
parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied Fsj outlier 
(BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in 
the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait 
associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and 
contemporary oxygen concentration) representing the strength of eutrophication. 

Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis 
methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in 
the more polluted lakes. 

Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in 
more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that 
the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of 
divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the 
number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of 
anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of 
heightened genomic divergence and therefore be classified as Est outlier loci. This study explores for the first time the 
potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence 
between coexisting species. 

Keywords: Coregonus, Population genomics. Adaptive radiation, Speciation reversal. Relaxed selection. Heterogeneous 
genomic divergence. Ecological speciation, Speciation-with-gene-flow, Eutrophication, Divergence hitchhiking 
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Background 

Understanding the evolutionary processes that control 
the diversification, coexistence and persistence of species 
has taken on a heightened urgency as human-induced 
environmental changes and stressors have increased in 
prevalence and the augurs of a new, anthropogenic mass 
extinction become more evident [1-3]. Changes in abi- 
otic and biotic environmental conditions will alter the 
direction, strength and form of selection that generates 
and maintains species differences and species coexist- 
ence, particularly among recently diverged species. How 
populations of such species respond to environmental 
change will depend on the severity, rate and duration of 
the change and the associated selection pressures experi- 
enced. Negative population growth is just one way spe- 
cies can go extinct; another scenario for species loss in 
the face of environmental change is speciation reversal 
[4]. Here, a reduction in environmental heterogeneity in- 
creases gene-flow and/or weakens disruptive selection, 
resulting in a break down in pre-zygotic and/or extrinsic 
post-zygotic reproductive isolation: leading to the fusion 
of previously distinct species through introgressive 
hybridization [5-8]. 

The effects of human-induced habitat alteration on 
the levels of intra- and interspecific diversity within 
adaptive radiations have generally been studied through 
the measurement of changes in the frequencies and dis- 
tributions of neutral genetic markers and/or putative 
adaptive trait values among and within the affected spe- 
cies. Examples from vertebrate radiations indicate a 
weakening of disruptive selection during the course of 
environmental perturbations and concomitantly an in- 
crease in hybridization between coexisting species, as 
the distribution of feeding resources and/or habitats for 
breeding becomes homogenized and prezygotic and ex- 
trinsic postzygotic reproductive isolation weakens. This 
process is accompanied by the complete or near 
complete fusion of previously discrete gene pools into 
single, albeit genetically diverse populations, and trait 
distributions losing multimodality [5,6,9]. 

During ecological speciation gene flow between diver- 
ging populations is expected to be restricted mainly at 
loci underpinning traits involved in divergent adaptation 
and areas of the genome in physical linkage with these 
loci; so called "genomic islands of divergence" [10,11]. 
Among recently diverged species, the balance between 
disruptive selection and homogenizing gene flow will de- 
termine the extent of genetic differentiation achieved 
and its genomic distribution [12,13]. This heterogeneous 
genomic divergence will be at its most extreme during 
moderate to high gene flow and strong exogenous dis- 
ruptive selection, predicted to characterize the early 
stages of sympatric/parapatric speciation scenarios 
[14,15]. As speciation progresses, genomic islands may 



increase both in number and extent, as inter-specific re- 
combination rates decrease, potentially enhancing re- 
productive isolation via divergence hitchhiking (DH) 
[10,13,16]. At later stages of ecological speciation, as many 
loci diverge among incipient species, genome-wide effect- 
ive gene flow will be reduced even at unlinked loci: diver- 
gence hitchhiking being superseded by genome hitchhiking 
(GH) [16]. Many studies have employed population gen- 
omic methods to identify genomic regions with unusually 
strong differentiation as holding potential candidate genes 
for disruptive selection between diverging taxa e,g, [17-21]. 
However little effort has been devoted to investigating how 
changes in natural selection regimes through anthropo- 
genic perturbation affect the occurrence and size of gen- 
omic islands of divergence between diverging (now 
converging) species, and nothing is published about this to 
our knowledge. 

The Alpine whitefish {Coregonus spp.) radiation is na- 
tive to around 40 lakes in the headwaters of three major 
European river drainages: the Rhine, Rhone and Danube 
[22]. Colonising following the glacial retreat (15,000 - 
20,000 years BP), this monophyletic group has diverged 
in-situ into at least seven species flocks, centred on sin- 
gle large lakes or previously connected sister lakes 
[23-25]. The constituent species of these lake or super- 
lake flocks are generally each other s closest relatives, in- 
dicating intralacustrine speciation [25]. Historically, 
Coregonus species richness varied from one to six en- 
demic species per lake. Coexisting species differ in traits 
linked to feeding and reproductive ecology, with very 
similar ecotypes having evolved repeatedly in several 
lakes [23,25,26]. 

During the 20th Century, many of the lakes in which 
the Alpine Coregonus radiation occurs have seen unpre- 
cedented anthropogenic impacts. Run-off from agricul- 
ture and insufficient sewage treatment massively 
increased the nutrient input into many lakes, leading to 
eutrophication of these naturally oligotrophic lakes. The 
increase in biologically available phosphorus had several 
effects on lake ecology, changing the community struc- 
ture of phytoplankton and their zooplankton grazers, 
therefore affecting the prey resource base for whitefish 
species [8,27]. Increased productivity also resulted in 
increased bacterial decomposition, reducing levels of 
available oxygen in the hypolimnion [28,29]. This de- 
oxygenation is most severe at the sediment-water inter- 
face, where whitefish egg development occurs [29]. A 
consequence of this is the reduction of the available 
depth gradient partitioned between whitefish species 
during spawning and concomitantly an increase in rates 
of hybridization between shallow and deep spawning 
species [8,30,31]. In more extreme cases de-oxygenation 
renders deeper areas of the lake anoxic, potentially redu- 
cing the available area and diversity of foraging habitat 
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[29]. During this period of unprecedented environmental 
change, the Alpine whitefish radiation experienced a loss 
of around 38% of its pre-eutrophication species richness 
[8]. This diversity loss appears to have been driven by a 
combination of negative population growth combined 
with the genetic collapse of coexisting species through 
speciation reversal Indeed by comparing historical 
population samples taken before large-scale eutrophica- 
tion with samples taken from the same species post- 
eutrophication, we see a clear signal of lower levels of 
neutral genetic and phenotypic differentiation between 
coexisting species within lakes that had high levels of eu- 
trophication, in comparison to species from lakes that 
have never been strongly eutrophic [8] . 

Here we present the first genome scan analysis that 
explores the impact of anthropogenic environmental 
change on heterogeneous genomic divergence in an 
adaptive radiation. We studied whitefish species from 
nine Alpine lakes, comprising five phylogenetically inde- 
pendent Coregonus adaptive radiations and one species 
assemblage arising from recent anthropogenic secondary 
contact. We employ genome scans of 183 individuals 
comprising 851 AFLP markers. We subject these AFLP 
loci to FsT outlier detection to identify candidate regions 
in the genome that likely evolved under disruptive selec- 
tion, and logistic regression analyses methods to link 
candidate regions to adaptive traits that they potentially 
influence. For all analysis methods, we consistently find 
a trend for more outlier and more trait-associated loci in 
those lakes that experienced lower levels of eutrophica- 
tion. Here we provide the first genetic evidence for the 
weakening of the genomic signature of disruptive selec- 
tion and divergent adaptation that are predicted to cause 
speciation reversal following human-induced habitat 
modification. 

Methods 

Study system & sampling 

We sampled Coregonus species from nine different lakes, 
each containing between two and five coexisting species. 
Individuals were collected between 2004 and 2006 on 
their respective spawning grounds using overnight gill- 
netting. Information about the species and lakes sampled 
can be found in Table 1 and Figure 1. Species were 
assigned to feeding ecotype using their mean gill-raker 
number, a correlate of feeding ecotype [32]: very low 
gill-raker numbers (VLGR <25 gill-rakers), low (LGR 
25-30), medium (MGR 31-35) and high (HGR >35). 
We used both the historical maximum total dissolved 
phosphorus concentration (Pmax (Ptot) ((^g 1^)) and con- 
temporary minimum oxygen concentration at depth 
(Min. O2 (mg 1'^)) as indicators of the intensity of an- 
thropogenic nutrient pollution: both variables potentially 
show different aspects of the speciation reversal process. 



relating to the duration and intensity of habitat disrup- 
tion and its effects on genomic divergence. P^ax was 
used as a surrogate for the maximum level of pollution 
experienced by each lake system and therefore to meas- 
ure the impact on interspecific genomic divergence 
when pollution was at its most severe. Min. O2 was used 
to gauge the extent to which lake systems have currently 
returned to their pre-eutrophication state, and therefore 
represents the potential ongoing effects of relaxed selec- 
tion on genomic divergence. Both P^ax and Min. O2 
have been shown to strongly correlate (positively and 
negatively respectively) with a loss in intralacustrine 
Coregonus species diversity [8]. 

Lakes Neuchatel and Biel are connected by ca, 8 Km 
of river and their three endemic Coregonus species form 
part of a monophyletic group with other nearby single- 
species lakes [23,25]. The five Coregonus species 
inhabiting connected {ca. 5 km of river) lakes Thun and 
Brienz, form a monophyletic clade too, except C fatioi, 
which appears to have had some genetic input from 
Lake Constance Coregonus species [23,25]. The endemic 
species flock of Lake Lucerne consists of potentially five 
coexisting species [Lundsgaard et al, unpublished obser- 
vations, 26]. For this study no samples of two recently 
discovered forms, intermediate in size and gill-raker 
counts, were available. Lake Constance consists of two 
major basins, Obersee (476 Km^) and Untersee (63 
Km^). Previously with up to six reported species, heavy 
eutrophication has lead to the extinction of at least one 
species, the VLGR C. gutturosus [8,22]. Four species co- 
exist in the Obersee, with the C. sp, "weissfelchen" found 
in the Untersee basin [23,25]. Lakes Walen and Zuerich 
are connected by ca, 17 Km of river, the two Coregonus 
species endemic to these lakes, together with whitefish 
from two nearby, single-species lakes, form an independ- 
ent, monophyletic radiation. Additionally from Walen and 
Zuerich, one or two poorly known species, may also have 
become extinct [8,26]. The two coexisting whitefish spe- 
cies of Lake Maggiore are introduced; C. sp, "lavarello" is 
most closely related to C. sp, "Zugerbalchen" from Lake 
Zug [23,25], while C. sp, "bondella", groups closest to C. 
zugensis from Lake Lucerne [25]. This species pair repre- 
sents a case of recent secondary contact, C. sp, "lavarello" 
and C. sp, "bondella" being introduced ca, 1860 and 1950 
respectively [33]. 

AFLP analysis 

Total genomic DNA was extracted using the Wizard 
Genomic DNA purification kit (Promega, Madison, WI). 
A genome scan of putatively independent loci was then 
carried out using AFLPs, following a protocol modified 
from [34], with combined restriction/ligation steps, 
using £coRI and Msel restriction enzymes. Selective 
PCRs were performed with 14 different primer-pairs 
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Table 1 Geographical location, lake surface area, lake trophic state, taxonomy, spawning/feeding ecotype designation 



and the number of individuals of each species included, for all species flocks analysed 

Lake Co-ordinates Surface Area (km^) Pmax (Ptot) (MQ I^) Min. O2 (mg 1'^) Species Ecotype # Individuals 

Neuchatel 46°53 N; 217.9 50 3.2 C. palaea LGR, W 7 

6°5 E C. condidus MGR, W 7 

14 

Biel 47°04N; 39.3 132 0.5 C. poloeo LGR, W 6 

7°1 E C. confusus MGR, W 7 

13 

Thun 46°41 N; 48.4 21 7.6 C. olpinus VLGR, S 7 

7°42 E C. sp. "balchen" LGR, W 6 

C. fatioi MGR, W 8 

C. sp. "felchen" MGR, W (Deep) 7 
C. a/6e//LiS HGR, S 7 

35 

Brienz 46°43 N; 29.8 17 9.3 C. sp. "balchen" LGR, W 4 

7°58 E C. sp. "felchen" MGR, W 7 

C. a/6e//us HGR, S 5 

16 

Lucerne 47°01 N; 114 34 3.4 C. sp. "bodenbalchen" LGR, W 10 

8°22 E C. zugensis HGR, S/W 9 

C. nobilis HGR, S 7 

26 

Constance 47°36 N; 539 87 6.0 C. sp. "weissfelchen" LGR, W 7 

9°26 E C. orenicolus LGR, W 6 

C. sp. "Alpenrhein" MGR, W 5 

C. macrophthalmus HGR, W 6 

C. wartmanni HGR, W (Pelagic) 6 

30 

Walen 47°07 N; 24.1 26 6.1 C. dup/ex LGR, W 8 

9°12E Cheglingus HGR, S/W 9 

17 

Zuerich 47°1 3 N; 88.4 119 1.1 C. duplex LGR, W 1 1 

8.43 E C. heglingus HGR, W 7 

18 

Maggiore 45°58 N; 212.5 37 n/a C. sp. "lavarello" MGR, W 8 

8°4 E C. sp. "bondella" HGR, W 6 

14 



Feeding ecotype: VLGR = Very low mean gill-raker number (<25 gill-rakers), LGR = low (25-30), MGR = medium (31-35) and HGR = high (>35). Reproductive 
ecotype: W = winter spawner (October-January), S = summer spawner (July-September). P^ax = peak historical levels of phosphorus enrichment. Min. 
O2 = contemporary minimum oxygen concentration at depth. 



(ACG-CAC, ACT-CAC, ACC-CAG, ACG-CAG, ACT- 
CAG, ACT-CAT, ACA-CGC, ACC-CGC, ACT-CGC, 
ACT-CGG, ATG-CTC, ATG-CTG, AAG-CTT, ACG- 
CTT). Fragment separation was carried out by elec- 
trophoresis on a CEQ-8000 automatic sequencer 
(Beckman-Coulter). AFLP traces were then screened and 



scored for peak presence/absence using GENEMARKER 
1.85 (SoftGenetics, LLC, State College, PA). Individual 
traces were automatically scored and then bins were 
visually inspected and manually adjusted. Only AFLP 
loci between 55-400 base pairs were scored. To check 
levels of marker reproducibility 30 individuals were 
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Figure 1 Geographic locations of the species floclcs studied within Switzerland. Map showing the number of species included, their 
sampling location and their representative phenotypes (species names listed from above to below). (A) Lake Neuchatel: (C. palaea, C. candidus), 
(B) Lake Biel: (C. poloeo, C. confusus), (C) Lake Thun: (C. olpinus, C. sp. "balchen", C. fotioi, C. sp. "felchen", C. olbellus), (D) Lake Brienz: (C. sp. 
"balchen", C. sp. "felchen", C. albellus), (E) Lake Lucerne: (C. sp. "bodenbalchen", C. zugensis, C. nobilis), (F) Lake Zuerich: (C. duplex, C. heglingus), (G) 
Lake Walen: (C. duplex, C. heglingus), (H) Lake Constance: (C. sp. "weissfelchen", C. arenicolus, C. sp. "Alpenrhein", C. macrophthalmus, C. wartmanni), 
(I) Lake Maggiore: (C. sp. "lavarello", C. sp. "bondella"). Yellow circles denote native species flocks. Red circle denotes anthropogenic origin and 
secondary contact. 



repeated (16% of the total analysed samples). The level 
of reproducibility per locus was then calculated for each 
of the 14 genotyped primer-pairs [35]. 

Genomic outlier analysis 

All outlier analyses and logistic regressions were run on 
the nine lake flock data sets separately. This avoided po- 
tential violations of the demographic models employed 
in the outlier detection methods that could be caused by 
hierarchical structuring of genetic variation [36,37], 
present in the Alpine whitefish radiation [25]. Addition- 
ally, with this analysis design, we were able to compare 
paired lake flock data sets with similar evolutionary his- 
tories and patterns of ecological diversification, but dif- 
fering in their nutrient enrichment history: Lake Walen 
(very low level of nutrient enrichment) vs. Lake Zuerich 
(high). Lake Biel (high) vs. Lake Neuchatel (moderate to 
high), and Lake Thun (low) vs. Lake Brienz (very low) 
(Figure 2). 

Two different Fsx-based approaches to search for se- 
lection candidate loci, suitable for dominant markers 
were used: DFDIST_D [38] and BAYESCAN [39]. 
DFDIST estimates the probability of a locus to be under 



selection according to its observed Fst and He, com- 
pared to simulated neutral distributions computed for 
an average neutral Fst observed in the empirical data set 
(the trimmed mean Fst)> effective population size and 
minimum levels of polymorphism [38]. In contrast, 
BAYESCAN directly estimates the probability of each 
locus being subject to selection by comparing the likeli- 
hoods of neutral and selective models [39]. As the power 
to detect stabilising selection will be low among recently 
diverged species, only candidate loci showing unusually 
high levels of genetic differentiation, implicating disrup- 
tive natural selection, were counted [38]. Fst outlier ana- 
lyses were first run on each multi-species lake data set 
(global comparisons among all species sampled in a 
lake), followed by all possible pairwise comparisons be- 
tween individual species for lakes with more than two 
species. 

In the DFDIST analysis the mean Fst, used for simulat- 
ing a null distribution of Fst values for loci evolving 
under neutrality, was estimated using a trimmed mean 
of 0.3 [38]. The null distribution was generated using 
50,000 simulated loci. In the empirical AFLP data set the 
maximum allele frequency was set to 0.98. AFLP loci 
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Figure 2 Yearly trends in total phosphorus concentration 


(ng r^) for sampled lakes of the Alpine Coregonus radiation. 


Phosphorus measurements were taken during lake overturning 


(modified from [8]). Colours of trend lines refer to lake identity: 


Neuchatel (pink), Biel (green), Thun (light blue), Brienz (red). Lucerne 


(dark blue), Constance (purple), Walen (grey), Zuerich (black) and 


Maggiore (orange). 



with FsT values above the upper 95th quantile of the null 
distribution were considered candidate loci for being 
under disruptive natural selection. 

In the BAYESCAN analysis, model parameter estimation 
was automatically tuned during the pilot runs (ten runs, 
length 2,000). A 10,000 iteration burn-in was found to be 
sufficient for convergence between MCMC chains. Sample 
sizes of 10,000 and a thinning ratio of 50 were used as sug- 
gested by [39]. Prior to the analysis all non-polymorphic 
loci (fixed presence or absence alleles) were removed as 
this could bias the analysis (M. Foil pers. comm.). Loci 
were considered as candidates for being under divergent 
selection if their logio Bayes Factor scores were > 0.5 ("Sub- 
stantial" on Jeffreys scale of evidence) [39]. 

To understand the effects of loci under selection on 
the phylogenetic relationships within the Alpine white- 
fish radiation, three different sets of Neighbour-joining 
(NJ) population-based trees were created: (1) using all 
genotyped AFLP loci (835 loci), (2) a subset comprising 
only outlier loci identified by DFDIST and BAYESCAN 
in either global or pairwise comparisons (96 loci), (3) a 
subset comprising only neutral loci remaining following 
removal of outlier loci (739 loci). Neis genetic distances 
between populations were calculated in AFLP-SURV 1.0 
[40]. NJ trees were reconstructed from these distances 
using the program NEIGHBOR in the software package 
PHYLIP 3.67, with 1,000 bootstrap replicates [41]. To 
quantify the levels of parallel ecotypic trends among lakes 
in the different loci subsets, following Allender et al [42], 
Mantel randomizations were performed for neutral and 
outlier loci subsets for a matrix of Neis genetic dis- 
tances and a design matrix encoding feeding ecotype. In 
the design matrix, species belonging to the same eco- 
type were encoded D = 0, D = 1 encoding assignment to 



different ecotypes. As the Maggiore species are non- 
native, they were excluded from this analysis. Mantel 
randomizations of 10,000 permutations were carried 
out in ARLEQUIN 3.5 [43]. 

To assess the statistical independence of individual 
outlier loci from each other, levels of linkage disequilib- 
rium (LD) were estimated for different AFLP loci sub- 
sets. For each species, genotyped AFLP loci were split 
into exclusive subsets: neutral and outlier. The outlier 
subsets comprised all loci identified as outliers by 
DFDIST analysis, during global-level comparisons within 
their respective species flock. LD (r^) values were calcu- 
lated in the program MULTILOCUS 1.3 [44]. Pairwise 
fd values were calculated for all possible pairs of loci 
within each neutral and outlier loci subset, for each spe- 
cies. The significance of differences in mean pairwise LD 
between neutral and outlier loci subsets, across popula- 
tions, was tested using paired ^-tests. 

Divergent trait measurement 

The first two trait categories measured, spawning depth 
and time, are related to reproductive ecology [26,31]. (1) 
Spawning date influences the time of emergence of the 
fry and therefore may influence survival through differ- 
ent mechanisms such as food availability, predation or 
over- wintering mortality [45-48] and has been show to 
transgressively segregate in whitefish hybrids [48,49]. 
Importantly, spawning date also influences temporal re- 
productive isolation between species [26]. The spawning 
date for each individual was taken as the day of capture 
on spawning site, adjusted so that the earliest caught fish 
(summer spawning) were caught on day 1, therefore 
those fish caught x days after the first fish were given 
the spawning date of Fish caught on the same date of 
different calendar years would be given the same 
spawning date. (2) Spawning depth is likely related to 
physiological or behavioural adaptations to depth [50], 
to different predation pressures [51], or to differences in 
egg size and development [48]. Also spawning depth in- 
fluences spatial reproductive isolation between species 
[8,26,31]. Spawning depth for each individual was esti- 
mated by finding the average depth over which the gill 
nets were set during each sampling event. 

Finally we measured three trait categories primarily 
related to feeding ecology: gill-raker number, linear mor- 
phometric distances and geometric shape analyses. (3) 
Gill-raker number has been shown to be correlated with 
feeding niche in Scandinavian whitefish [52] and to have 
a high heritability [32]. Gill-raker numbers were counted 
on the first gill arch from the left side of each fish. (4-5) 
Overall morphological variation among fish was mea- 
sured using (4) linear morphological distances (LM) and 
(5) geometric morphometries of body shape (GM). 
Shape aspects are potentially targets of disruptive 
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selection, as revealed in Qst-Fst comparisons in lakes 
of the Alpine radiation [31] and show significant 
phenotype-environment correlations in Scandinavian 
Coregonus species flocks [53]. (4) Fifteen different LM 
measurements, based on those described in Chouinard, 
Pigeon & Bernatchez [54], were taken using digital calli- 
pers on 1090 individuals, comprising 38 whitefish popu- 
lations/species from 22 lakes in the Alpine region: to 
capture the major axes of morphological divergence 
across the radiation, (Figure 3b-d). Size correction of the 
measured traits was performed using common principal 
components analysis (CPCA), followed by Burnaby's 
back projection method (BBPM) through the CPCBP 
package [55] in R (R Development Core Team, 2004). 
Principal component analysis (PCA) was then performed 
on the transformed data in SPSS 16 (SPSS Inc., Chicago 
Illinois, USA). The resulting PC 1-3 scores of the AFLP- 
genotyped individuals from the nine lake flocks were 
taken for logistic regression analysis. For PC 1 (30% of 
variation) traits with positive loadings were mainly body 
length traits (AFL, PFL, DFL, HL), whereas head mea- 
sures particularly jaw width traits (LJW, UJW, UJL, ED, 
SEL) had negative loadings. PC 2 (24%) had positive 
loadings for head measures (HL, ED, SEL, EH, LJL, UJL) 
and negative loadings for traits: AFL, PFL, DFL, TD, 
LJW, UJW. PC 3 (10%) had strongly positive loadings 
for eye shape (ED, EH) and negative loadings for: SEL, 
UJL, INB, SNW. (5) GM landmark-based analysis of 
body shape variation was carried out on digital images 
of 1545 fish comprising 41 populations/species from 22 
lakes, using 13 homologous landmarks (Figure 3a). 
Digital landmark placement was carried out in TPSDIG2 
(http://life.bio.sunysb.edu/morph/). All further shape 
analysis was carried out in MORPHOJ 1.02d [56]. Gen- 
eralized least-squares Procrustes super-imposition was 
carried out on landmark configurations. AUometric size 
correction was then carried out by the pooling of 
species-specific regressions of Procrustes coordinates 
against individual standard length. Following PCA, PC 1 
was found to be mostly artifactual (associated with 
dorso-ventral flexion of dead fish) and was discarded. 
Therefore individual PC scores for PC 2 (13%), PC 3 
(9%) and PC4 (9%) were used for further analysis. Shape 
changes along PC 2 (13%) were associated with changes 
in head length. PC 3 (9%) was associated mainly with 
the amount of up/down rotation of the head, and 
changes in body depth. PC 4 (9%) was mainly associated 
with changes in the length of the caudal peduncle and 
body depth. 

Association between genomic outlier loci and putative 
adaptive traits 

Univariate logistic regressions were carried out in the 
program MATSAM 2Beta [57]. This approach analyses 




UJW 



1.LJW 2. IND 



Figure 3 Landmarks and distance-based measurements used in 
the analysis of Coregonus body shape. 13 landmarks used to 
analyse body shape with geometric morphometries (a). Fifteen linear 
morphometric measurements taken for each individual: (DFL) length 
to dorsal fin, (PFL) length to pelvic fin, (AFL) length to anal fin, (TD) 
tail depth, (HL) head length, (ED) eye diameter, (SNL) snout length, 
(EOL) eye to operculum length, (EH) eye height, (LJL) lower jaw 
length, (UJL) upper jaw length, (INB) intra-orbital distance, (IND) 
intra-nostril distance, (LJW) lower jaw width, (UJW) upper jaw width 
(b-d). Modified from [34]. 



the relationship between alleles at a given AFLP locus 
and specific values of an explanatory variable, such as 
putative adaptive traits [57]. For each AFLP locus, the 
genotype of individuals was regressed against the indi- 
vidual values for each of the five putative adaptive trait 
variables measured: spawning date, spawning depth, gill- 
raker number, PC score on principal axes of morpho- 
logical variation resulting from LM measurements and 
GM analysis of shape. Each suite of genotype/phenotype 
regressions were calculated separately for each of the 
nine species flock data sets. Individuals with missing 
phenotypic data were removed from the logistic regres- 
sion for that specific trait. The adaptive trait values for 
each individual, used in logistic regressions, can be 
found in Additional file 1: Table SI. The significance of 
the association between the frequency of an allele at a 
locus and the value of each trait was assessed using the 
G-statistic at the 0.01, 0.05 and 0.1 significance levels, 
with Bonferroni correction using the number of AFLP 
loci used in each regression analysis. 
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Tests of association between heterogeneous genomic 
divergence and eutrophication 

To test if the distribution over lakes of the detected gen- 
omic signal of disruptive selection (number of outlier 
loci and trait-associated loci) was robust to analysis 
method (DFDIST, BAYESCAN, MATSAM), correlations 
between the numbers of loci detected by each analysis 
method were carried out. To test the hypothesis that eu- 
trophication causes a decrease in the number of identi- 
fied outlier/trait-associated loci, one-tailed regressions 
were calculated of the number of loci detected in global 
comparisons within each species flock, against max- 
imum historical phosphorus concentration (Pmax) and 
current minimum oxygen concentration at depth (Min. 
O2) of the lake the species flock inhabits. This test was 
repeated for each of the three methods employed 
(DFDIST, BAYESCAN, MATSAM). 

As a comparison the effect of eutrophication on neu- 
tral markers was also tested by linear regression of P^ax 
and Min. O2 against Est values from each lake flock 
calculated from three different marker sets: (1) 
microsatellites (maximum within-lake pairwise Est 
value), (2) all AELP loci (highest within-lake weighted 
mean Est in DEDIST from pairwise comparisons), (3) 
neutral AELP loci (highest within-lake weighted mean 
Est in DEDIST from pairwise comparisons, following 
the removal of all candidate loci identified by the 
DEDIST, BAYESCAN, MATSAM analyses). Lake flock 
microsatellite Est values resulting from between ten 
(Maggiore only) or eleven unlinked loci were taken from 
[8,58], Hudson et aL unpublished observations. We also 
calculated pairwise correlations between the three inde- 
pendent estimates of Est to see how congruent their sig- 
nals were. All statistical analyses were carried out in 
SPSS 16 (SPSS Inc., Chicago Illinois, USA). 

Results 

The total number of AELP loci scored per lake flock var- 
ied from between 355-405 (DEDIST) and 216-306 
(BAYESCAN) (Table 2). The number of whitefish geno- 
typed per flock ranged from 13 to 35 (Tables 1, 2), de- 
pending largely on the number of constituent species. 
Eor species-pairwise analyses the number of loci geno- 
typed ranged from between 334-380 (DEDIST) and 
177-247 (BAYESCAN), and the number of individuals 
from 9-19. Average per locus levels of reproducibility 
across the fourteen primer pairs exceeded 95%. Identity 
codes, the total number of loci (both polymorphic and 
fixed) and average reproducibility per primer pair can be 
found in Additional file 2: Table S2. 

FsT outlier scans 

The DEDIST analyses detected disruptive selection can- 
didate loci in every species flocks. The largest number of 



outliers was found in the Lake Thun species assemblage 
(22) (Table 2, Additional file 3: Table S3a). Most loci 
detected (64 of 81) fell within the 95-99% quantiles. 
Loci with Est values above the 99% quantile were found 
in L. Thun (four), L. Zuerich (three), L. Brienz (two), L. 
Lucerne (two), L. Walen (two), L. Biel (one) and L. Con- 
stance (one). Of the four loci found to be outliers in 
multiple lakes, three loci showed parallel allelic trends 
among similar ecotypes. CAGC-235 was an outlier in 
the Thun, Lucerne, Zuerich and Maggiore species flocks. 
Eor this locus prevalence of the presence allele was re- 
duced in Zuerich LGR C. duplex, Lucerne LGR C. sp, 
"bodenbalchen", Thun MGR C. fatioi and Maggiore 
MGR C. sp, "lavarello". The presence allele of CGAG- 
186 was at a higher prevalence in Brienz MGR C. sp, 
'felchen" and Thun/Brienz HGR C. alhellus. At TGTG- 
183 there was a decreased prevalence of the presence al- 
lele in LGR C. sp. 'balchen" from both lakes Brienz and 
Thun and also in Thun MGR C, fatioi, Pairwise DEDIST 
analyses (Additional file 3: Table S3b) revealed 36 add- 
itional loci not recovered in the global analyses. Six of 
these loci were identified as outliers in several pairwise 
species comparisons from different lakes, with three 
showing parallel trends among similar ecotypes. Locus 
CAGC-342 had a high prevalence of the presence allele 
in HGR L. Brienz C. albellus and HGR L. Constance C. 
wartmanni, TGTG-271 showed parallel trends in HGR 
Lucerne C. zugensis and HGR Maggiore C. sp. 
"bondella". Finally the presence allele of CGAG-355 had 
a reduced frequency in MGR Constance C. sp. 
"Alpenrhein" and MGR Maggiore C. sp. "lavarello". 

From the global lake-level BAYESCAN analyses, we find 
disruptive selection candidate loci between the constituent 
species in every lake except Biel and Constance (Table 2, 
Additional file 4: Table S4a). 15 out of the 20 outlier 
loci identified had Bayes factor scores of "substantial" 
on Jeffreys scale of evidence (> 0.5 - 1). One locus scor- 
ing "strong" (> 1 - 1.5) was found in each of Thun, Brienz 
and Zuerich flocks and two in Maggiore. Similar to the 
DFDIST analyses, CAGC-235, was detected as an outlier in 
multiple lakes. The presence allele for this locus had a 
low prevalence in the Thun MGR C. fatioi and Lucerne 
LGR C. sp. 'bodenbalchen". Pairwise BAYESCAN analyses 
(Additional file 4: Table S4b) revealed five additional outlier 
loci not recovered in the global analyses at lake level. None 
of these loci were outliers in more than one species pair. 

Overall the DFDIST analysis identified over four times 
more outlier loci than the BAYESCAN analysis (81 vs. 
20). Similarity of identity of candidate loci identified with 
either method (global Est analysis at lake level) was low- 
est in lakes with low average mean Est (for Lakes Con- 
stance and Biel zero loci were found in the BAYESCAN 
analyses, but four and two loci respectively in DFDIST 
analyses) and/or with more than two species coexisting 
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Table 2 Number of selection candidate loci identified using DFDIST and BAYESCAN within individual lake flocks 
(global comparisons) 



Species Flock 


# Species 


# Individuals 


# Poly. Loci 


# Candidate loci 
DFDIST 


Total 


# Candidate loci 
BAYESCAN 


Total 


Shared Loci * 


Yo Similarity 










> 0.95 > 0.99 




> 0.5 


> 1 








Neuchatel 


2 


14 


241 


5 


0 


5 


2 


0 


2 


2 


40 


Biel 


2 


13 


262 


1 


1 


2 


0 


0 


0 


0 


0 


Thun 


5 


35 


306 


18 


4 


22 


3 


1 


4 


4 


18 


Brienz 


3 


16 


216 


9 


2 


11 


1 


1 


2 


2 


18 


Lucerne 


3 


26 


266 


9 


2 


11 


3 


0 


3 


2 


18 


Constance 


5 


30 


279 


4 


1 


5 


0 


0 


0 


0 


0 


Walen 


2 


17 


229 


2 


2 


4 


2 


0 


2 


2 


50 


Zuerich 


2 


18 


230 


5 


3 


8 


3 


1 


4 


4 


50 


Maggiore 


2 


14 


223 


11 


2 


13 


1 


2 


3 


3 


23 








Total 


64 


17 


81 


15 


5 


20 







Presented are the numbers of species, individuals and polymorphic AFLP loci comprising each species flock data set. Identified DFDIST outlier loci are grouped 
according to how many fell between the 95% and 99% quantiles (> 0.95), and how many fell above the 99% quantile (> 0.99). BAYESCAN outlier loci are grouped 
by their Bayesian posterior probability of being under selection (Logio Bayes factor scores): "substantial" 0.5 to 1 (> 0.5), "strong" 1 to 1.5 (> 1). Also shown: the 
number of AFLP loci identified as selection candidates by both outlier detection methods and the overall similarity between sets of candidate loci recovered in 
DFDIST and in BAYESCAN, for each species flock. 



(L. Brienz, L. Thun and L. Lucerne; only 18% of loci 
identified by both detection methods) (Table 2). Devia- 
tions were in all cases due to DFDIST detecting higher 
numbers of Fst outliers. 

Consensus trees from the three different subsets of 
AFLP loci showed high levels of congruence (Additional 
file 5: Figure SI). Compared to the all loci consensus tree, 
the outlier loci subset tree showed little obvious increase 
in ecotype clustering. Design matrices of feeding ecotype 
designation explained slightly more variance in genetic 
distances between species when the latter were calculated 
from the outlier loci subset rather than the neutral loci 
subset (r = 0.07 and r = 0.04 respectively), but none of 
these Mantel randomizations were significant (0.004%, 
P = 0.12 and 0.001%, P = 0.22 respectively). 

The mean values of pairwise LD for neutral and outlier 
loci classes, across species, were not significantly differ- 
ent (mean difference = 0.026, t23 = L591, p = 0.125). 
Overall outlier loci subsets generally had lower mean 
pairwise LD within species than neutral loci subsets 
(0.077 V5. 0.105). The range in mean pairwise LD values 
across species was greater for the outlier loci class (0 - 
0.333 V5. 0.077 - 0.157), however variance around 
population means were qualitatively similar for both out- 
lier and neutral loci classes (Additional file 6: Table S5, 
Additional file 7: Figure S2). 

Association between genomic outlier loci and putative 
adaptive traits 

Across all lake flocks, a total of 43 significant associa- 
tions between putative adaptive traits and AFLP markers 
were found (Table 3, Additional file 8: Table S6). The 



divergent trait with the largest numbers of significant as- 
sociations with AFLP alleles, across all lake flocks, was 
gill-raker number (19 loci), followed by LM PC 1 (jaw 
widths) (eight), spawning depth (six), spawning date 
(four), GM PC 2 (head length), LM PCs 2 (eye and head 
size) and 3 (skull widths) (each two). No significant asso- 
ciations were found for GM PCs 3 (head position) and 4 
(caudal peduncle length). Within lake flocks the largest 
number of significant associations were found in Lake 
Thun (17), followed by L. Lucerne (nine), L. Walen (six), 
L. Brienz (five), L. Zuerich (four), L. Constance (two) 
and none in lakes Neuchatel, Biel and Maggiore. 

Four loci showed significant parallel trends of associ- 
ation with adaptive traits in multiple lake flocks. As in 
the candidate loci screen, CAGC-235 again showed par- 
allel allelic trends; the presence allele was significantly 
associated in Lakes Lucerne and Brienz with lower gill- 
raker numbers. The CTAC-99 presence allele was sig- 
nificantly associated with increased gill- raker numbers in 
lakes Brienz and Thun. The presence allele of CGAG- 
186 was found in individuals with high gill-raker num- 
bers in L. Brienz and individuals with more negative LM 
PC 3 scores (reduced skull widths) in L. Thun. The pres- 
ence allele of C TAG -173 had a significant association 
with earlier spawning dates in L. Thun and with lower 
gill-raker numbers in L. Lucerne. 

The total number of trait-associated loci found using 
MATSAM, also detected as Fst outlier loci was higher 
for DFDIST (26) than for BAYESCAN analyses (11). 
However the average percentage of loci that were both 
outlier and trait-associated, were very similar between 
the sets of loci identified with DFDIST (22%), and 
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Table 3 Number of AFLP loci revealing a significant association with each putative adaptive trait, for each species flock 

# Significant Loci 



Species Flock 


GR 


LM PC 1 


LM PC 2 


LM PC 3 


GM PC 2 


Date 


Depth 


Total 


Neuchatel 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0 


Biel 


0/ 0/ 0 


n/a 


n/a 


n/a 


0/ 0/ 0 


0/ 0/ 0 


0/ 0/ 0 


0 


Thun 


2/4/4 


3/ 5/6 


0/ 0/0 


1/2/2 


0/ 0/ 1 


0/0/3 


0/ 0/ 1 


17 


Brienz 


0/ 2/4 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/ 1 


0/ 0/0 


5 


Lucerne 


1/2/5 


0/ 0/0 


0/ 0/ 1 


0/ 0/0 


0/ 0/0 


0/ 0/0 


2/ 3/3 


9 


Constance 


0/ 0/ 1 


0/ 0/ 1 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


2 


Walen 


0/ 1/ 1 


0/ 1/ 1 


0/ 0/ 1 


0/ 0/0 


0/ 0/ 1 


0/ 0/0 


0/ 0/2 


6 


Zuerich 


0/4/4 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


4 


Maggiore 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0/ 0/0 


0 


Total 


19 


8 


2 


2 


2 


4 


6 


43 



Numbers in adaptive trait columns refer to the total number of loci found at each level of significance using the G-statistic {p - < 0.01/ p - < 0.05/ p - < 0.1). 
Adaptive traits used: gill-raker number (GR), linear morphometric PGA scores (LM PG 1, PG 2, PG 3), geometric morphometric PGA scores (GM PG 2), spawning date 
and spawning depth, n/a = analyses that could not be carried out due to missing data. GM PG 3, PG 4 are not represented as no significant associations were 
found for these variables. 



BAYESCAN (21%) (Additional file 9: Table S7). Overall 
the lake flock with the highest number of selection can- 
didate loci significantly associated with adaptive trait 
values was Thun (12 loci). Lakes Brienz, Lucerne and 
Zuerich came next highest, each having 4 loci, with 
Walen having two such loci, the remaining species flocks 
having none. 

Tests of association between heterogeneous genomic 
divergence and eutrophication 

The number of selection candidate and trait-associated 
loci identified at the species flock level through the three 
different analyses methods were significantly positively 
correlated in two of the three comparisons: DFDIST ver- 
sus BAYESCAN (rg = 0.717, P < 0.05) and DFDIST ver- 
sus MATSAM (rg = 0.763, P < 0.05). Results from 
BAYESCAN and MATSAM also tended to be positively 
correlated (rg = 0.570, P = 0.1). All linear regressions of 
the effects of Pmax on the number of selection candidate 
loci and significant trait-locus associations showed nega- 
tive, albeit marginally non-significant relationships (DFDIST 
global comparisons ^ = 0.315, P = 0.058; BAYESCAN glo- 
bal comparisons ^ =0.155, P = 0.148; MATSAM ^ = 0.239, 
P = 0.09; all P values one-tailed: Figure 4). The regression 
effect of Min. O2 on the number of selection candidate loci 
and trait-locus associations was positive, but weaker than 
that of Pmax and non-significant (DFDIST global compari- 
sons = 0.263, P = 0.097; BAYESCAN global compari- 
sons =0.009, P = 0.413; MATSAM / = 0.236, P = 0.111; 
Additional file 10: Figure S3). 

Tests of association between neutral genomic divergence 
and eutrophication 

In Linear regressions, P^ax values were strongly nega- 
tively correlated with current lake flock-wide global Est 



for neutral microsatellites, all AFLP markers and the 
neutral AFLP subsets (r^ = 0.38-0.65, P = 0.008-0.079; 
Table 4, Figure 5). Levels of Min. O2 were significantly 
positively correlated with lake flock-wide global Est and 
explained even more variation than P^ax (^^ = 0.56-0.76, 
P = 0.005-0.033; Additional file 11: Figure S4). Esti- 
mates of Est from all three datasets were highly posi- 
tively correlated across lakes (r = 0.906-0.975, d.f. = 8, 
P <0.001). 

Discussion 

During early stages of ecological speciation, gene flow 
between incipient species in sympatry/parapatry is 
expected to be restricted mainly at loci underlying adap- 
tive traits and associated physically linked loci: resulting 
in increased levels of genetic differentiation in these re- 
gions of the genome, relative to unlinked loci evolving 
under neutrality. In this study we find candidates for loci 
evolving under disruptive natural selection in each of 
the examined lake species flocks of the Alpine Coregonus 
adaptive radiation. We also find significant associations 
between many of these candidate loci and putative adap- 
tive traits, supporting the action of disruptive natural se- 
lection in the origin and coexistence of species within 
this radiation. Some of the lakes have recently under- 
gone anthropogenic homogenization of whitefish habitat 
through eutrophication, and all analysis methods show 
negative trends in the number of disruptive selection 
candidate loci (DFDIST, BAYESCAN) and the number 
of significant loci-adaptive trait associations (MATSAM) 
found among coexisting species, in relation to levels of 
historical eutrophication within the lakes the whitefish 
species flocks inhabit. This is congruent with patterns of 
genetic differentiation shown by neutral microsatellite 
and AFLP loci in the same species flocks, indicating the 



Hudson et al. BMC Evolutionary Biology 2013, 13:108 
http://www.bionnedcentral.conn/1471 -21 48/1 3/1 08 



Page 11 of 1 7 



a 



I 

1 10- 

S 



t 

u 

v 

c 

1 = 







R» = 0.315 




#TlnMi 






• Magglore 






#^ "■■#l,ucerne 


#Zuerlch 




9 Neuchatel Constance 
Walen 


Blel 


0 


25 50 75 100 


125 




Maximum Phosphorus Cone, (pg l-*) (P^ax) 






R' = 0.239 




•Thun 






% Lucerne 






^ Walen 


^Zuerich 




• 

Magglore Neuchatel Constance 


Biel 
•••^ 


0 


25 50 75 100 
Maximum Phosphorus Cone, (pg 1 ^) (P, 


125 

ax) 







R»= 0.155 


Thun 




#Zuerlch 


Lucerne^ 

Magglore 






Brienz 

• • • 

Walen Neuchatel 








Constance 


Blel 



Maximum Phosphorus Cone, (pg l-*) (Ph^x) 



Figure 4 The effects of historical phosphorus concentration (Pmax) on the number of candidate loci identified. Linear regressions of tine 
effects of liistorical pliospliorus concentration (Pmax) on tlie number of candidate loci identified witliin lal<e flocl<s in (a) DFDIST analyses, (b) 
BAYESCAN analyses and (c) the number of significant trait-loci associations in MATSAM. 



weakening of reproductive isolation during eutrophica- 
tion. We hypothesize that the pattern of reduced numbers 
of candidate loci found in species flocks that have experi- 
enced more severe anthropogenic eutrophication, could be 
due to two distinct but non-exclusive genomic processes 
brought about by increased interspecific hybridization and 
a concomitant increase in interspecific recombination. 
In the first, increased interspecific recombination 
erodes genomic islands of elevated genetic differenti- 
ation surrounding the target loci of disruptive natural 
selection, akin to a reversal of divergence hitchhiking, 
reducing the likelihood that AFLP restriction sites will 
fall within these regions. In the second, large chromo- 
somal segments, including intact genomic islands of 
divergence, pass between the genomes of hybridizing 
species, equalizing allele frequencies and reducing 
heightened Fst levels associated with genomic islands 
and therefore reducing the power to detect outlier 
candidate loci within those segments. Our results pre- 
sented here provide the first genomic evidence that eu- 
trophication may lead to the relaxation of divergent 
selection that in turn permits increased levels of inter- 
specific, introgressive hybridization, driving speciation 
reversal. 



Evidence for the role of disruptive natural selection in 
shaping intralacustrine diversity 

Previous phylogenetic and population genetic analyses 
have shown that the Alpine Coregonus radiation is 
monophyletic and that its constituent lake flocks form a 
nested series of independent radiations that arose rap- 
idly, in-sitUy within the last 15,000 years [25]. During the 
course of an adaptive radiation, disruptive natural and 
sexual selection regimes experienced by species in con- 
trasting environments drive adaptive divergence and 
speciation [59]. Given the co-occurrence of multiple, 
ecologically divergent species within the Alpine lake 
flocks and the overall youth of the radiation, heteroge- 
neous genomic divergence between divergent species is 
expected [14]. Natural selection against hybrids with 
trait values or trait combinations maladaptive in either 
parental environment will maintain stronger differenti- 
ation in the face of gene flow in genomic islands of ele- 
vated interspecific genetic differentiation surrounding 
loci under disruptive selection than in regions invisible 
to selection [12]. We find candidate loci displaying un- 
usually strong patterns of genetic differentiation between 
paired coexisting species in each of the species flock 
studied. 
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Table 4 Maximum within-lake Fst values among 
sympatric species derived from different genetic markers 

Species Flock AFLP Microsatellite 



All Neutral 



Neuchatel 


0.077 


0.028 


0.05* 


Biel 


0.073 


0.035 


0.024t 


Thun 


0.153 


0.107 


0.252t 


Brienz 


0.229 


0.185 


0.269t 


Lucerne 


0.123 


0.074 


0.13* 


Constance 


0.087 


0.039 


0.071* 


Walen 


0.095 


0.049 


0.13* 


Zuerich 


0.095 


0.045 


0.044* 


Maggiore 


0.173 


0.097 


0.182+ 



Comparison of tlie maximum witliin-lal<e pairwise Fst values for eacli lake 
flock. All = Fst value derived using all genotyped AFLP loci. Neutral = Fst value 
derived following removal of all AFLP selection candidate loci. Data taken 
from: ^ Vonlanthen et al. (2012), t Bittner (2009), t Hudson et al. 
(unpublished observations). 



Given the recent, common origin of the Alpine radi- 
ation, parallel trends in candidate loci between similar 
ecotypes in different lake flocks could be expected. 
Shared standing genetic variation and similar levels of 
genetic co-variation among traits increases the likelihood 
that the same loci and alleles are involved in the evolu- 
tion of equivalent phenotypes in different lake flocks 
[60]. We find support for this notion, albeit the number 
of loci that are repeatedly detected as outliers and/or as- 
sociated with putative adaptive traits is moderate. Also 
there is no obvious increase in the phylogenetic group- 
ing of species by ecotype in trees constructed solely from 
identified outlier loci, and matrices coding for ecotype 
designation in Mantel tests explained non-significant 
levels of variation in interspecific genetic distances gen- 
erated solely from identified outlier loci: though stiU 
more than the variation explained in tests using neutral 
AFLP loci. The moderate number of parallel allelic 
trends likely reflects that many of the candidate loci we 
identified, are linked with varying degrees of physical 
proximity to loci under disruptive natural selection, ra- 
ther than being the target loci of disruptive selection 
themselves [13,61]. During speciation within the inde- 
pendent whitefish species flocks, the action of recombin- 
ation, eroding blocks of linkage disequilibrium along the 
genome, and the time since origination of the beneficial 
alleles, could mean that different genetic backgrounds 
containing different hitchhiker AFLP outlier loci are 
driven to, or near to fixation in similar ecotypes, among 
species flocks: despite the true target locus of selection, 
that hitchhiking loci are physically linked too, being 
identical in each paralleled speciation event [60]. There- 
fore parallel trends in candidate loci between similar 
ecotypes would likely be restricted to hitchhiker loci very 



tightly linked to the actual loci under selection and/or 
parallel population divergences of recent, shared evolu- 
tionary origins. 

From our study, locus CAGC-235 showed parallel al- 
lelic trends in more than two lake systems. The presence 
allele for this locus has a reduced prevalence in lower 
gill-raker ecotypes in lakes Maggiore, Thun (C. fatioi 
only), Zuerich, Lucerne and Brienz (C. sp, "balchen" 
only). Therefore there is the potential that the alternative 
allele for CAGC-235 is tightly linked to the actual target 
locus for natural selection. Five of the remaining loci 
support the role of recombination in breaking down the 
physical linkage between homologous candidate loci and 
the target of selection over evolutionary time: the paral- 
lel trends in these loci being shared only between the 
most recently separated lake flocks. Three loci were 
found in both lakes Thun and Brienz (CGAG-186, 
TGTG-183, CTAC-99), which form part of the same 
monophyletic super-lake radiation and given the geo- 
graphical proximity of the lakes, were likely the last of 
the studied species flocks to separate. CGAG-355 and 
TGTG-271 show paraUel allelic trends between intro- 
duced L. Maggiore species and their likely source popu- 
lations in either lakes Lucerne and Constance [25]. 

Anthropogenic perturbations and heterogeneous 
genomic divergence 

The genomic extent of levels of increased genetic differen- 
tiation around the target loci for disruptive natural selec- 
tion will depend on the balance between the strength of 
selection and gene flow [12,13]. Two distinct, but to some 
degree complimentary, mechanisms have been proposed 
to describe how genetic differentiation spreads through 
the genomes of species diverging under speciation-with 
-gene-flow scenarios. Under divergence hitchhiking (DH), 
strong disruptive selection acting on a locus reduces inter- 
specific recombination over relatively restricted areas of 
the genome in physical linkage with the selective locus. In 
this way beneficial alleles under weaker selection can be 
recruited for adaptive divergence through their "uplift" 
onto genomic islands of divergence surrounding more 
strongly selected loci. This reduces the chances of their 
loss to interspecific recombination and facilitates the 
growth of genomic islands of divergence [13,62]. Under 
genome hitchhiking (GH), disruptive selection reduces the 
average effective migration across the entire genome, 
thereby facilitating differentiation throughout the genome 
irrespective of levels of physical linkage to selective loci 
[63,64]. Feder, Egan & Nosil [64] describe a four-phase 
model of speciation-with-gene-flow. In this model, DH is 
most important in facilitating genomic differentiation at 
early to intermediate stages of speciation (phase 2). At 
later stages, as more loci are recruited by disruptive nat- 
ural selection, GH will become the dominant mechanism 
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driving genomic differentiation (phase 3). The relatively 
recent origin of the species flocks of the Alpine Coregonus 
radiation (<15,000 years) suggests that DH could play an 
important role in the divergence seen between species. 
Under DH, genomic islands of divergence around selective 
loci are expected to increase in size as speciation pro- 
gresses. This in turn would increase the likelihood that in- 
dividual AFLP loci will fall within these areas of higher 
genetic differentiation and be classed as outlier candidate 
loci. Therefore as genomic divergence progresses under 
DH, the proportion of outlier candidate loci amongst a 
random genomic sample of AFLP loci should increase: 
relative to a similar genome scan conducted at an earlier 
time point in the speciation process of the same species 
pair [12,13]. This prediction depends on the existence of 
relatively extensive genomic regions of low interspecific 
recombination around selected loci. In dwarf/normal spe- 
cies pairs in the closely related Coregonus clupeaformis 
species complex, effective migration was found to be 
reduced for relatively large distances around selective loci 
(~ 5 cM) and that the extent of phenotypic divergence be- 
tween coexisting species pairs was positively correlated 
with levels of linkage disequilibrium around these loci 
[48,65]. However, in the species flocks of the Alpine white- 
fish radiation more genomic research will be needed to 



characterize the actual target loci of selection and define 
the effective limits of surrounding islands of divergence, 
taking into account factors which could potentially en- 
hance {e,g, epistasis) or constrain {e,g, selection from 
standing genetic variation) the size of these genomic areas: 
before extensive DH can be confirmed [13,63]. 

In the lakes of the Alpine Coregonus radiation, previ- 
ous research has shown that historical levels of an- 
thropogenic nutrient enrichment negatively correlated 
with contemporary interspecific genetic differentiation at 
neutral markers, and positively correlated with the over- 
all loss of species and functional phenotypic diversity. 
Direct comparisons of neutral genetic differentiation 
from samples taken from the same species flocks, pre- 
and post- (large-scale) eutrophication, corroborate this 
pattern, with strong declines found in a highly 
eutrophied lake relative to more lightly polluted lakes, 
despite similar levels of genetic differentiation occurring 
in these species flocks prior to the peak of eutrophica- 
tion [8]. Our present study corroborates these results, 
with genetic differentiation at selectively neutral subsets 
of AFLP loci also being strongly negatively correlated 
with levels of anthropogenic nutrient enrichment. Focus- 
ing on areas of the genome evolving under the influence 
of disruptive natural selection, we hypothesized that the 
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numbers of candidate loci identified within a lake flock 
would be also negatively correlated with each lake s his- 
torical maximum level of eutrophication. As the lake 
habitats become increasingly homogeneous during eu- 
trophication, disruptive selection between habitats would 
wane, reducing the efficacy of barriers to interspecific 
gene flow and resulting in increased interspecific recom- 
bination that would, in turn, either reduce the size of 
genomic islands of divergence, reversing divergence 
hitchhiking or result in the transfer of entire chromo- 
somal segments containing intact genomic islands be- 
tween the genomes of hybridizing species, equalizing 
allele frequencies at potential Fst outlier loci. However, 
whilst we find that numbers of candidate loci in lakes 
with less strong eutrophication histories were indeed 
high (lakes Thun, Brienz, Lucerne) and were lowest in 
the most eutrophic lake (Biel), the relationship with 
phosphorus was marginally non-significant. This was 
due to certain lakes having fewer candidate loci 
(Neuchatel, Walen) and others lakes having more 
(Zuerich), than would be predicted given their past levels 
of eutrophication. 

Variation in the levels of ecological divergence be- 
tween coexisting species, among species flocks, prior to 
environmental perturbations could potentially affect the 
relative numbers of candidate loci found among lake 
flocks today. If ecological divergence and therefore de- 
velopment of RI was impeded by high levels of gene 
flow, weaker disruptive selection and/or temporal fluctu- 
ations in selection regimes, this would affect the number 
and extent of genomic islands of divergence that occur 
between coexisting species naturally and how susceptible 
they are to environmental disturbance when it occurs 
and subsequently the extent of speciation reversal, rela- 
tive to similarly aged species flocks from more ecologic- 
ally divergent/stable lake environments [4,66]. The two 
pairs of coexisting species from recently connected lakes: 
Walen and Zuerich, and Biel and Neuchatel, represent 
two distinct super-lake clades within the overall Alpine 
radiation [25]. Lakes Zuerich and Biel have a similar 
history of high eutrophication (Pniax=119 (ig/l and 
132 (ig respectively), L. Neuchatel moderate (50 (ig/l) 
and L. Walen low (26 (ig/l). Using gill raker number as a 
surrogate for ecological diversity among coexisting spe- 
cies, L. Zuerich had a greater pre-eutrophication differ- 
ence in average gill-raker number between species (26.7 
& 35.2) than lakes Biel (27.4 & 33.2) or Neuchatel (27.5 
& 32.2), similar to L. Walen (25.8 & 35.4) [8,29]. Follow- 
ing heavy eutrophication, despite evidence of increased 
neutral gene flow and low microsatellite Est values simi- 
lar to Biel and Neuchatel [8], a higher number of candi- 
date loci was found and relatively high phenotypic 
differentiation has been maintained between L. Zuerich 
species, despite much more severe levels of pollution 



relative to Neuchatel. Additional support for this 
hypothesis may come from the closely related. North 
American C. clupeaformis -complex. The coexisting 
Coregonus species of Indian pond show high levels of 
phenotypic differentiation and the largest number of 
AFLP candidate loci yet found in the C. clupeaformis ra- 
diation, yet have relatively low neutral genetic differenti- 
ation [18,67]. Intriguingly, Indian Pond has phosphorus 
concentrations three times higher than other studied 
North American Coregonus lakes and the highest post- 
summer decreases in oxygen with depth, suggesting this 
lake and its endemic Coregonus species may have experi- 
enced eutrophication in its recent past [68]. 

In pristine environments, as the number of species 
within a radiation increases so should the number of se- 
lection candidate loci and trait associated loci identified, 
as the number of axes of divergent adaptation among 
coexisting species increases [65]. In the Lake Thun spe- 
cies flock, five species coexist and the highest numbers 
of candidate loci/trait associated loci found in this study 
were identified within this species flock. In contrast, the 
similarly species -rich, yet heavily polluted L. Constance 
species flock had some of the lowest numbers of can- 
didate loci/trait associated loci identified in the Alpine 
radiation. L. Constance previously held up to six 
coexisting species, but following eutrophication, the 
VLGR ecotype, C. gutturosus went extinct and among 
the remaining species a large decrease in interspecific 
neutral genetic differentiation occurred [8]. This may il- 
lustrate a potential trend in speciation reversal, in that 
young lake radiations with higher species diversity are 
intrinsically more vulnerable, in response to a given level 
of environmental change, to speciation reversal than 
less species-rich systems. Within evolutionarily young, 
species-rich radiations, ecological and reproductive 
niche breadths may be narrower and more tightly 
packed compared to related, less diverse radiations, per- 
haps through partial subdivision of the ancestral niche: 
potentially making them more sensitive to relaxation of 
disruptive natural selection regimes. As environmental 
factors imparting RI begin to break down the increased 
potential directions of gene flow between multiple 
coexisting species could in turn exacerbate speciation re- 
versal. This hypothesis could help to explain rapid col- 
lapses in diversity within other young radiations such as 
the Lake Victoria cichlids [5] and Great Lake ciscoes [9]. 

The strategy we have adopted for the present study 
was to investigate all possible Alpine whitefish species 
flocks, including most of the endemic species present: 
each species flock comprising a distinct data point along 
a spectrum of anthropogenic habitat perturbation from 
low to severe. Overall, the limited number of individuals 
per species included in our analyses certainly places lim- 
itations on the power to identif)^, among the set of loci 
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sampled, all those that might be evolving under disrup- 
tive selection between coexisting species. Also, whilst 
mean values of linkage disequilibrium were not higher 
among identified outlier loci than among neutral loci, 
suggesting that detected outliers have a rather dispersed 
distribution in the genomes of Alpine whitefish species, 
the true number of separate linkage groups that our can- 
didate loci comprise remains unknown. Additionally, we 
cannot rule out entirely that another, unidentified eco- 
logical factor, covarying with the extent of eutrophica- 
tion, may be a stronger driver of the pattern seen in the 
number of outlier loci identified amongst Swiss whitefish 
species flocks. In this regard, our study stands as a ne- 
cessary first step, allowing future studies to home in on 
the actual target loci of disruptive selection and their 
function, to uncover the genomic architecture of the 
phenotypic traits they affect and to measure the fitness 
effects of alleles at these loci within the alternate niches 
of the different species, and how eutrophication may 
affect these. Given the number and diversity of species 
flocks covered by our sampling scheme, which explicitly 
takes into account the hierarchical structuring of genetic 
variation within the Alpine whitefish radiation; we have 
the power to look at the relative numbers of candidate 
loci identified and formulate hypotheses as to the factors 
influencing variation in this number. Therefore despite 
the fact that while some of the individual candidate loci 
identified in this study, may have elevated genetic differ- 
entiation through processes other than contemporary 
disruptive natural selection, such as historical selective 
sweeps [69], demography [70] or their location in 
chromosomal regions with reduced recombination [21]; 
the recent, common origin of the species flocks of the 
Alpine whitefish radiation means that these potential 
biases are very unlikely to lead to the trends in candidate 
loci we find across lake flock data sets. Furthermore by 
using two different Fst outlier detection methods, along- 
side logistic regressions, combining population-based 
and individual-based approaches respectively, not only 
can we identify adaptive traits important in niche differ- 
entiation and the loci potentially underpinning them, 
but also that the high congruence in the results pro- 
duced by these different methodologies suggests that the 
overall patterns in candidate loci revealed by our study 
are real. 

Conclusions 

Within each of the postglacial species flocks of the Al- 
pine whitefish radiation we find evidence for the action 
of disruptive natural selection in the shaping of within- 
lake species diversity through the identification of Est 
outlier loci and/or significant adaptive trait/candidate 
loci associations. We also find, similar to the extent 
of interspecific genetic differentiation at neutral loci, a 



negative correlation between the number of identified 
FsT-outlier loci identified and significant adaptive trait/ 
candidate loci associations, and historical levels of an- 
thropogenic nutrient enrichment within each lake. We 
hypothesize that the habitat disruption caused by eu- 
trophication weakens disruptive natural selection re- 
gimes resulting in increased gene flow between species. 
The resulting interspecific recombination either breaks 
down genomic islands of heightened genetic differenti- 
ation surrounding the target loci of disruptive natural se- 
lection or entire chromosomal segments, including 
intact genomic islands, pass between species' genomes, 
homogenizing allele frequencies. Our study presents the 
first genomic evidence that speciation reversal is associ- 
ated with the relaxation of disruptive selection: contrib- 
uting to our understanding of the mechanisms that 
control this important driver of biodiversity loss. 

Additional files 



Additional file 1: Table SI. Individual phenotypic trait values. 

Additional file 2: Table S2. AFLP primer pair information. 

Additional file 3: Table S3. Identity of AFLP loci classified as Fst 
outliers following DFDIST analyses. 

Additional file 4: Table S4. Identity of AFLP loci classified as Fsj 
outliers following BAYESCAN analyses. 

Additional file 5: Figure SI. Population-based NJ consensus trees 
created using different subsets of AFLP loci. 

Additional file 6: Table S5. Mean pairwise LD values from neutral and 
outlier loci classes for each whitefish species. 

Additional file 7: Figure S2. Distribution of pairwise linkage 
disequilibrium values for neutral and outlier loci classes within each 
whitefish species. 

Additional file 8: Table S6. Identity of AFLP loci with significant 
associations with putative adaptive phenotypic traits identified by 
MATSAM. 

Additional file 9: Table S7. Congruence between Global dataset 
evidence for disruptive selection and association with adaptive traits 
within species flocks. 

Additional file 10: Figure S3. Linear regressions of the effects of 
minimum lake oxygen concentration at depth (Min. O2) on the number 
of candidate loci identified within lake flocks in (a) Dfdist analyses, 
(b) BayeScan analyses and (c) the number of significant trait-loci 
associations in MatSAM. 

Additional file 11: Figure S4. Linear regressions of the effects of 
minimum oxygen concentration at maximum depth (Min. O2) on genetic 
differentiation in different sets of loci: (a) neutral microsatellite loci, (b) all 
AFLP loci (c) 'neutral' AFLP loci. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

AGH wrote the paper, contributed to conception and design of the study, 
sampled fish, generated adaptive trait and genetic data, and carried out 
most of the subsequent analyses. PV sampled fish, generated environmental 
lake data and contributed to analyses and writing. EB contributed to 
analyses and writing. OS conceived and designed the project, supervised the 
project, and contributed to data analyses and writing. All authors read and 
approved the final manuscript. 



Hudson et al. BMC Evolutionary Biology 2013, 13:108 
http://www.bionnedcentral.conn/1471 -21 48/1 3/1 08 



Page 16 of 17 



Acknowledgements 

Many thanks to R. Mueller, D. Bittner, B. Lundsgaard-Hansen, M. Kugler, J. 
Muggli, B. Poll!, P. Volta, D. Gerdaux and A. Champigneulle for assisting in 
sampling. Also thanks to K. Lucek, S. Mwaiko, A. Sivasundar, V. Schneider, D. 
Joyce, C. Klingenberg and B. M. Bolker for technical advice and assistance. 
Constructive suggestions from D. Roy, M. Barluenga, I. Magalhaes and P. 
Nosil greatly improved this manuscript. This study was funded by the 
University of Bern and the EAWAG Handlungsfeld project: Aquadiverse. 

Author details 

^Division of Aquatic Ecology & Macroevolution, Institute of Ecology & 
Evolution, University of Bern, Baltzerstrasse 6, CH 3012 Bern, Switzerland. 
^Department of Fish Ecology & Evolution, EAWAG Centre of Ecology, 
Evolution and Biogeochemistry, Seestrasse 79, CH 6047 Kastanienbaum, 
Switzerland. ^Department of Ecology & Environmental Science, Umea 
University, SE 901 87 Umea, Sweden. ^Department of Biology, Reed College, 
3203 SE Woodstock Blvd, Portland, OR 97202, USA. 

Received: 19 March 2013 Accepted: 20 May 2013 
Published: 30 May 2013 

References 

1. Smith TB, Bernatchez L: Evolutionary change in human-altered 
environments. Mol Ecol 2008, 17:1-8. 

2. Butchart SHM, Walpole M, Collen B, van Strien A, Scharlemann JPW, Almond 
REA, Baillie JEM, Bombard B, Brown C, Bruno J, Carpenter KE, Carr GM, 
Chanson J, Chenery AM, Csirke J, Davidson NC, Dentener F, Foster M, Galli 

A, Galloway JN, Genovesi P, Gregory RD, Hockings M, Kapos V, Lamarque J-F 
, Leverington F, Loh J, McGeoch MA, McRae L, Minasyan A, Morcillo MH, 
Oldfield TEE, Pauly D, Quader S, Revenga C, Sauer JR, Skolnik B, Spear D, 
Stanwell-Smith D, Stuart SN, Symes A, Terney M, Tyrrell TD, Vie J-C, Watson 
R: Global biodiversity: indicators of recent declines. Science 2010, 
328:1164-1168. 

3. Pereira HM, Leadley PW, Proenga V, Alkemade R, Scharlemann JPW, Fernandez- 
Manjarres JF, Araujo MB, Balvanera P, Biggs R, Cheung WWL, Chini L, Cooper 
HD, Gilman EL, Guenette S, Hurtt GC, Huntington HP, Mace GM, Oberdorff T 
Revenga C, Rodrigues P, Scholes RJ, Sumaila UR, Walpole M: Scenarios for 
global biodiversity in the 21st century. Science 2010, 330:1496-1501. 

4. Seehausen 0: Conservation: Losing biodiversity by reverse speciation. 
CurrBiol 2006, 16:334-337. 

5. Seehausen 0, van Alphen JJM, Witte F: Cichlid fish diversity threatened by 
eutrophication that curbs sexual selection. Science 1997, 277:1808-181 1. 

6. Taylor EB, Boughman JW, Groenenboom M, Sniatynski M, Schluter D, Gow 
JL: Speciation in reverse: morphological and genetic evidence of the 
collapse of a three-spined stickleback {Gasterosteus oculeotus) species 
pair. Mol Ecol 2006, 15:343-355. 

7. Hendry AP, Farrugia TJ, Kinnison MT: Human influences on rates of 
phenotypic change in wild animal populations. Mol Ecol 2008, 17:20-29. 

8. Vonlanthen P, Bittner D, Hudson AG, Young KA, Muller R, Lundsgaard- 
Hansen B, Roy D, Di Piazza S, Largiader CR, Seehausen 0: Eutrophication 
causes speciation reversal in whitefish adaptive radiations. Nature 2012, 
482:357-362. 

9. Todd TN, Stedman RN: Hybridization of ciscoes {Coregonus spp.) in Lake 
Huron. Can J Zool 1989, 67:1679-1685. 

1 0. Wu CI: The genie view of the process of speciation. J Evol Biol 2001 , 

14:851-865. 

11. Turner TL, Hahn MW, Nuzdin SV: Genomic islands of speciation in 
Anopheles gambiae. PLoS Biol 2005, 3:e285. 

12. Nosil P, Funk DJ, Ortiz-Barrientos D: Divergent selection and 
heterogeneous genomic divergence. Mol Ecol 2009, 18:375-402. 

13. Via S, West J: The genetic mosaic suggests a new role for hitchhiking in 
ecological speciation. Mol Ecol 2008, 17:4334-4345. 

14. Coyne JA, Orr HA: Speciation. Sunderland, MA: Sinauer Associates, Inc.; 2004. 

15. Mallet J, Meyer A, Nosil P, Feder JL: Space, sympatry and speciation. J Evol 
Biol 2009, 22:232-2341. 

16. Feder JL, Gejji R, Yeaman S, Nosil P: Establishment of new mutations 
under divergence hitchhiking and genome hitchhiking. Phil. Trans. R. Soc. 

B. 2012, 367:461-474. 

17. Wilding CS, Butlin RK, Grahame J: Differential gene exchange between 
parapatric morphs of Littorina saxatilis detected using AFLP markers. 
J Evol Biol 2001, 14:611-619. 



18. Campbell D, Bernatchez L: Generic scan using AFLP markers as a means 
to means to assess the role of directional selection in the divergence of 
sympatric whitefish ecotypes. Mol Biol Evol 2004, 21:945-956. 

19. Bonin A, Taberlet P, Miaud C, Pompanon F: Explorative genome scan to 
detect candidate loci fo adaptation along a gradient of altitude in the 
common frog (Rana temporaria). Mol Biol Evol 2006, 23:773-783. 

20. Nosil P, Egan SP, Funk DJ: Heterogenous genomic differentiation between 
walking-stick ecotypes: 'isolation by adaptation' and multiple roles for 
divergent selection. Evolution 2008, 62:316-336. 

21. Roesti M, Hendry AP, Salzburger W, Berner D: Genome divergence during 
evolutionary diversification as revealed in replicate lake-stream 
stickleback population pairs. Mol Ecol 2012, 21:2852-2862. 

22. Kottelat M, Freyhof J: Handbook of European freshwater fishes. Cornol, 
Switzerland: Publications Kottelat; 2007. 

23. Douglas MR, Brunner PC, Bernatchez L: Do assemblages of Coregonus 
(Teleostei: Salmoniformes) in the Central Alpine region of Europe 
represent species flocks? Mol Ecol 1 999, 8:589-603. 

24. Hudson AG, Vonlanthen P, Muller R, Seehausen 0: Review: The geography 
of speciation and adaptive radiation in Coregonines. Adv Limnol 2007, 
60:111-146. 

25. Hudson AG, Vonlanthen P, Seehausen 0: Rapid parallel adaptive radiations 
from a single hybridogenic ancestral population. Proc R Soc S 201 1, 

278:58-66. 

26. Steinmann P: Monographie der Schweizerischen Koregonen. Beitrag zum 
Problem der Entstehung neuer Arten. Spezieller Teil. Schweiz Z Hydrol 

1950, 12:340-491. 

27. Burgi HR, Heller C, Gaebel S, Mookerji N, Ward JV: Strength of coupling 
between phyto- and zoplankton in Lake Lucerne (Switzerland) during 
phosphorus abatement subseuent to a weak eutrophication. J Plankton 
Res 1999, 21:485-507. 

28. Gruber N, Wehrii B, Wuest A: The role of biogeochemical cycling for the 
formation and preseration of varved sediments in Soppensee 
(Switzerland). J Paleolinnnol 2000, 24:277-291 . 

29. Muller R, Stadelmann P: Fish habitat requirements as the basis for 
rehabilitation of eutrophic lakes by oxygenation. Fisheries Manag. Ecol. 
2004, 11:251-260. 

30. Numann W: The Bodensee: effects of exploitation and eutrophication on 
the Salmonid community. J. Fish. Res. Board Can. 1972, 29:833-847. 

31. Vonlanthen P, Roy D, Hudson AG, Largiader CR, Bittner D, Seehausen 0: 
Divergence along a steep ecological gradient in lake whitefish 
(Coregonus sp.). J Evol Biol 2009, 22:498-514. 

32. Bernatchez L: Ecological Theory of Adaptive Radiation. An Empirical 
Assessment from Coregonine Fishes (Salmoniformes). In Evolution 
lllunninated. Edited by Hendry AP, Stearns SC. Oxford: Oxford University 
Press; 2004:175-207. 

33. Berg A, Grimaldi E: Ecological relationships between planktophagic fish 
species in the Lago Maggiore. Verh. Int Ver Theor Angew. Limnol. 1966, 
16:1065-1073. 

34. Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T Homes M, Frijters A, 
Pot J, Peleman J, Kuiper M, Zabeau M: AFLP: a new technique for DNA 
fingerprinting. Nucleic Acids Res 1995, 23:4407-4414. 

35. Bonin A, Bellemain E, Eidesen PB, Pompanon F, Brochmann C, Taberlet P: 
How to track and assess genotyping errors in population genetics 
studies. Mol Ecol 2004, 13:3261-3273. 

36. Thornton KR, Jensen JD: Controlling the false-positive rate in multilocus 
genome scans for selection. Genetics 2007, 175:737-750. 

37. Excoffier L, Hofer T, Foil M: Detecting loci under selection in a 
hierarchically structured popultion. Heredity 2009, 103:285-298. 

38. Beaumont MA, Balding DJ: Identifying adaptive genetic divergence 
among populations from enome scans. Mol Ecol 2004, 13:969-980. 

39. Foil M, Gaggiotti 0: A Genome-Scan Method to Identify Selected Loci 
Appropriate for Both Dominant and Codominant Markers: A Bayesian 
Perspective. Genetics 2008, 180:977-993. 

40. Vekemans X: AFLP-SURV v. 1.0. Laboratoire de Genetique et Ecologie Vegetale, 
Universite Libre de. Distributed by the author: Bruxelles, Belgium; 2002. 

41. Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). 
Cladistics 1989, 5:164-166. 

42. Allender CJ, Seehausen 0, Knight ME, Turner GF, Maclean N: Divergent 
selection during speciaton of Lake Malawi cichlid fishes inferred from 
parallel radiations in nuptial coloration. P. Natl. Acad. Sci. USA 2003, 
100:14074-14079. 



Hudson et al. BMC Evolutionary Biology 2013, 13:108 
http://www.bionnedcentral.conn/1471 -21 48/1 3/1 08 



Page 1 7 of 1 7 



43. 



44. 



45. 



46. 



47. 



50. 



51. 



52. 



53. 



54. 



55. 



56. 



57. 



58. 



59. 



60. 



62. 



63. 



64. 



65. 



66. 



67. 



Excoffier L, Lischer HEL: Arelquin suite ver 3.5: A new series of programs 
to perform populaion genetic analyses under Linux and Windows. Mol 

Ecol Resour 2010, 10:564-567. 

Agapow P-M, Burt A: Indices of multilocus linkage disequilibrium. Mol. 
Ecol. Notes. 2001, 1:101-102. 

Gushing DH: Plankton production and year-class strength in fish 
populations - an update of the match mismatch hypothesis. Adv Mar 

Biol 1990, 26:249-293. 

Svensson E, Sinervo B: Experimental excursions on adaptive landscapes: 

Density-dependent selection on egg size. Evolution 2000, 54:1396-1403. 

Sogard SM: Size-selective mortality in the juvenile stage of teleost fishes: 

A review. B. Mar. Sci. 1 997, 60:1 1 29-1 1 57. 

Woods P, Muller R, Seehausen 0: Intergenomic epistasis causes 

asynchronous hatch times in whitefish hybrids, but only when parental 

ecotypes differ. J Evol Biol 2009, 22:2305-2319. 

Rogers SM, Bernatchez L: The genetic architecture of ecological 

speciation and the association with signatures of selection in natural 

lake whitefish (Coregonus sp. Salmonidae) species pairs. Mol Biol Evol 

2007, 24:1423-1438. 

Rogers SM, Gagnon V, Bernatchez L: Genetically based phenotype- 
environment association for swimming behavior in lake whitefish 
ecotypes (Coregonus clupeaformis Mitchill). Evolution 2002, 56:2322-2329. 
Weaver MJ, Magnuson JJ, Glayton MK: Habitat heterogeneity and fish 
community structure: Inferences from north temperate lakes. Am. Fish S. 
S. 1996, 16:335-346. 

Kahilainen K, Ostbye K: Morphological differentiation and resource 
polymorphism in three sympatric whitefish Coregonus lavaretus (L.) 
forms in a subarctic lake. J Fish Biol 2006, 68:63-79. 
Harrod G, Mallela J, Kahilainen KK: Phenotype-environment correlations in 
a putative whitefsh adaptive radiation, j /In/m Ecol 2010, 79:1057-1068. 
Ghouinard A, Pigeon D, Bernatchez L: Lack of specialization in trophic 
morphology between genetically differentiated dwarf and normal forms 
of lake whitefish (Coregonus clupeaformis Mitchill) in Lac de I'Est. 
Quebec. Can. J. Zool. 1996, 74:1989-1998. 

McGoy MW, Boll<er BM, Osenberg GW, Miner BG, Vonesh JR: Size correction: 
comparing morphological traits among populations and environments. 

Oecologia 2006, 148:547-554. 

Klingenberg GP: MorphoJ: an integrated software package for geometric 

morphometries. Mol. Ecol. Res. 201 1, 11:353-357. 

Joost S, Bonin A, Bruford MW, Despres L, Gonord G, Erhardt G, Taberlet P: A 
spatial analyss method (SAM) to detect candidate loci for selection: 
towards a landscape genomes approach to adaptation. Mol Ecol 2007, 
16:3955-3969. 

Bittner D: Gonad deformations in whitefish (Coregonus spp.) from Lake Thun, 
Switzerland - A population genetic and transcriptomic approach, PhD Thesis. 
Switzerland: University of Bern; 2009. 

Schluter D: The ecology of Adaptive Radiation. New York: Oxford University 
Press; 2000. 

Barrett RD, Schluter D: Adaptation from standing genetic variation. Trends 
Ecol Evol 2008, 23:38-44. 

Butlin RK: Population genomics and speciation. Genetica 2010, 

138:409-418. 

Via S, Gonte G, Mason-Foley G, Mills K: Localizing Fsj outliers on a QTL 
map reveals evidence for large scale genomic regions of reduced gene 
exchange during speciation-with-gene-flow. Mol Ecol 2012, 21:5546-5560. 
Feder JL, Gejji R, Yeaman S, Nosil P: Establishment of new mutations 
under divergence and genome hitchhiking. Philos Trans. R Soc. B: Biol Sci. 

2012, 367:461-474. 

Feder JL, Egan SP, Nosil P: The genomics of speciation-with-gene-flow. 

Trends Genet 2012, 28:303-306. 

Gagnaire P-A, Pavey SA, Normandeau E, Bernatchez B: The genetic 
architecture of reproductive isolation during speciation-with-gene-flow 
in lake whitefish species pairs assessed by RAD-sequencing. Evolution 

2013. doi:1 0.1 111/evo.l 2075. 

Nosil P, Harmon L, Seehausen 0: Ecological explanations for (incomplete) 

speciation. Trends Ecol Evol 2009, 24:145-156. 

Lu G, Bernatchez L: Correlated trophic specialization and genetic 

divergence in sympatric lake whitefish ecotypes (Coregonus 

clupeaformis): support for the ecological speciation hypothesis. Evolution 

1999,53:1491-1505. 



69. 



70. 



Landry L, Vincent WF, Bernatchez L: Parallel evolution of lake whitefish 
dwarf ecotypes in association with limnological features of the adpative 
landscape. J Evol Biol 2007, 20:971-984. 

Bierne N: The distinctive footprints of local hitchhiking in a varied 
environment and global hitchhiking in a subdivided population. 

Evolution 2010, 64:3254-3272. 

Gharlesworth B, Gharlesworth D, Barton NH: The effects of genetic and 
geographic structure on neutral variation. Annu Rev Ecol Syst 2003, 
34:99-125. 



doi:1 0.1 1 86/1 471 -21 48-1 3-1 08 

Cite this article as: Hudson et al.: Genomic signatures of relaxed 
disruptive selection associated with speciation reversal in whitefish. 

BMC Evolutionary Biology 201 3 1 3:1 08. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



